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p<| ' ABSTRACT 

o\ '. 

, Aims. We present a study of the highly ionized gas (plasma) associated with damped Lyman-a (DLA) systems at 2=2.1-3.1. 
i-H ' Methods. We search for Ovi absorption and corresponding Siiv, Civ, and Nv in a Very Large Telescope/Ultraviolet- Visible 
<^> Echelle Spectrograph (VLT/UVES) sample of 35 DLA systems with data covering Ovi at S/N>10. We then use optical depth 
C ■ profile comparisons and ionization modelling to investigate the properties, phase structure, and origin of the plasma. 

' Results. We report twelve DLAs (nine intervening and three at <5 000 kms" 1 from the QSO redshift) with detections of Ovi 
^ (—j absorption. There are no clear Ovi non-detections, so the incidence of Ovi in DLAs is between 34% (12/35) and 100%. Among 
O ,' these 12 DLAs, Civ and Siiv are seen whenever data is available, and Nv is detected in 3 cases. Analysis of the line widths 
I together with photoionization modelling suggests that two phases of DLA plasma exist: a hot, collisionally ionized phase (seen 
in broad O vi components) , and a warm, photoionized phase (seen just in narrow C IV and Si IV components) . The presence of 
inflows and/or outflows is indicated by individual Ovi and Civ components displaced from the neutral gas (either blueshifted 
| or redshifted) by up to 400 kms -1 . We find tentative evidence (98% confidence) for correlations between the DLA metallicity 
(measured in the neutral gas) and high-ion column density, and between the DLA metallicity and high-ion line width, as would 
be expected if supernova-driven galactic outflows rather than accretion produced the high ions. Using conservative ionization 
corrections, we find lower limits to the total hydrogen column densities in the hot (O vi-bearing) and warm (C iv-bearing) phases 

> ^ ■ in the range logTV^^' > 19.5 to > 21.1, and logA^ 1 ™ > 19.4 to > 20.9. On average, the hot and warm phases thus contain 

' ^40% and J>20% of the baryonic mass of the neutral phase in DLAs, respectively. 
Conclusions. If the temperature in the Ovi phase is « 10 6 K and so fom = O vi/O <C 0.2, the plasma can make a significant 
contribution to the metal budget at high redshift. Additional searches for Ovi in Lyman Limit Systems (QSO absorbers with 
17.0 < A?hi < 20.3) will be necessary to determine the total quantity of baryons and metals hidden in hot halos at z ~ 2. 
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1. Introduction 

Damped Lyman-a (DLA) systems, defined as those QSO 
absorbers with Hi column densities 2Vhi > 2 x 10 20 cm~ 2 , 
represent the largest reservoirs of neutral gas in the 
redshift range 0-5, and are believed to be the precur- 
sors i_to_j]radej^n_dj^ recent review 
bv IWolfe. Gawiser. fc Prochaskal [2005h. Detailed abun- 
dance studies of t h ese sy st ems at various redshifts 
dProchaska fe Wolfe] 2002 : Dessauges-Zavadskv et al. 
2004 120061) have provided a means to trace the process 



* Based on observations taken with the Ultraviolet and 
Visible Echelle Spectrograph (UVES) on the Very Large 
Telescope (VLT) Unit 2 (Kueyen) at Paranal, Chile, operated 
by ESO. 



of cosmic metal enrichment over a large fraction of the 
age of the Universe. 

In addition to the neutral gas phase, spectroscopic 
observations of DLAs have revealed both cold molecu- 
lar (T <, 300K ; IPetitiean et al.ll2000 HLedoux et al~l l2003: 
Srianand et alj 120051) and warm ionized (T ~ 10 4 K 



Lu et al l 11999: iLedoux et all Il998l : IWolfe fc Prochaska 



2000a) components, with a multi-phase structure resem- 



bling that seen in the Galactic interstellar medium (ISM). 
However, very little is known about the presence or prop- 
erties of a hot (T > 10 5 KjH ionized medium. Given the ev- 
idence for star formation in DLAs (see lWolfe et al.ll2005l ). 



1 There is some disagreement in the ISM and IGM literature 
about the meaning of the terms "warm" and "hot" . X-ray as- 
tronomers do not consider any gas below 10 6 K hot, and IGM 
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Type II supernovae should create observable regions of 
hot, shock- heated interstellar plasma. The separate pro- 
cess of accretion and shock-heating of intergalactic gas 
may also lead to the production of a hot ionized medium. 

The best ultra-violet (UV) lines available for study- 
ing hot interstellar plasma are the O VI AA1031, 1037 and 
Nv AA1238,1242 doublets. Ovi (Nv) peaks in abun- 
dance at T = 3(2) x 10 5 K under collisional ioniza- 



tion e^[uili^riu : m_J^IF^ conditions ([Sutherland fe Dopita 
19931; iGnat fe Sternberg||2006l ). Surveys of O VI absorption 
have been used to establish the presence of a hot ionized 
medium in the Milky Way and its surrounding network 
of high- velocity clouds (Savage et al.ll2003l : ISembach et al 
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Fox. Savage. fe Wakkerl l2006h . Additionally, Ovi 



can be formed by photoionization, as is reported for nar 
row O VI absorbers in the intergalactic medium (IGM) 
at z_ £ 2 dCarswell. Schave. fe Kim 2002t Bergeron et al 



20021: Bergeron fc Herbert-Ford 120051 : (Levshakov et al 
20031 ). Distinguishing between photoionized and colli- 



sionally ionized O vi absorbers is therefore important in 
searching for a hot ionized medium in DLAs. 

O vi absorption has been found before in Lyman Limit 
Systems (QSO absorbers with 17.0 < logA^i < 20.3; 
Kirkman fc Tvtledll997l Il999h but no systematic search 
for Ovi absorption has yet been conducted in DLAs. In 
this paper we present such a search. In Sect. 2 we discuss 
the data acquisition and reduction, O vi identification, and 
absorption line measurements. In Sect. 3 we present the 
spectra and compare the line profiles, and in Sect. 4 we 
discuss ionization processes. In Sect. 5 and 6 we discuss 
the total plasma content in DLAs and its contribution to 
the cosmic density 51. In Sect. 7 we briefly speculate on 
the origins of the O vi phase. We summarize our study in 
Sect. 8. 



2. Data acquisition and handling 

2.1. Observations 

The VLT/UVES DLA sample currently consists of 123 
DLA (logA^Hi > 20.3) and sub-DLA (19.7 < logA^i < 
20.3) systems. All data were acquired with the Ultraviolet- 
Visible Echelle Spectrograph (UVES) at the 8.2 m Very 
Large Telescope Unit 2 (VLT/UT2) at Paranal, Chile in 
the years 2000 to 2006. The data reduction was performed 
using the interactive pipeline described in iBallester et all 
(l2000h: full details o f the reduction procedures are given in 



Ledoux et al.l (|2003l ). The rebinned pixel size is «2 km s 1 
and the data have a spectral resolution (FWHM) of 
6.6 kms- 1 (i?=45000). 



astronomers now have adopted the term "warm-hot" to de- 
scribe the regime with 10 5 to 10 7 K. However, throughout this 
paper we use the ISM traditions of "warm" to imply T ~ 10 4 K 
and "hot" to mean T > 10 5 K. 



2.2. Identifying vi absorption in DLAs 

O vi systems are only accessible from the ground at 
z <; 2, where the transitions become redshifted enough 
to pass the atmospheric cutoff near 3000 A. They lie in 
the Lyman-a forest, the series of intervening Hi absorp- 
tion line s whose number density increases rapidly above 
z = 2.5 (|Kim et al. 1997 ). There is therefore a high level 
of confusion in separating O VI from H I interlopers, that 
becomes more difficult with increasing redshift. We thus 
adopted a series of systematic steps to identify our sample 
of DLA Ovi absorbers, as follows. 

(1) We selected from the UVES DLA sample the 35 
DLAs with data covering the O vi wavelength range at 
S/N>10 per resolution element^. 

(2) We immediately rejected 19 cases where the O vi 
lines arc clearly blended, e.g. by continuous absorption 
troughs extending over several hundred kms -1 . 

(3) We searched for Ovi components over a range of 
±500 kms -1 relative to the redshift of the neutral gas in 
the DLA. If absorption lines were found where the doublet 
ratio t(0 vi A1031) /r(0 vi A1037)=2 in all pixels through 
the line profile, we made a preliminary identification of 
O vi. Alternatively, if one O vi line was blended but the 
other O vi line showed a candidate absorption line, we 
looked at the Civ data (which exists in most cases), and 
if the candidate O vi absorber showed a similar line profile 
to the C iv absorber, we treated the O vi as a preliminary 
identification. 

(4) We investigated whether the preliminary O vi 
identifications could be caused by intervening Lyman-a, 
Lyman-/? or Lyman-7 forest absorbers, by looking in each 
case for corresponding absorption in the other Lyman 
series lines and in Civ. If no blending lines were iden- 
tified, the preliminary Ovi detections became real and 
entered our sample. Four preliminary O vi identifications 
were identified as being due to intervening Hi lines. We 
are confident that no Lyman-/? or Lyman-7 forest in- 
terlopers remain in our sample, since in these cases the 
corresponding Lyman-a would be seen in the spectrum. 
It is more difficult to decontaminate Lyman-a absorbers 
(which are more numerous) since the corresponding higher 
order Lyman lines often fall outside the observed wave- 
length range. Therefore, it is possible that some contami- 
nation by Lyman-a forest absorbers remains in our data. 

The problem of reliably discriminating between true 
high-rcdshift Ovi absorbers and intervening Lyma n for- 
est systems is disc ussed in iBurles fc Tvtler ( 1996 ) and 
ISimcoe et al. ( 2002 ). Our situation is favorable in that we 
look for O vi at the redshift z a b s where an absorption line 
system is already known to exist, not in a blind search 
over redshift. We can crudely assess the probability that a 
given O VI component is caused by blending using the ob- 
se rved density of Ly man-q forest lines. The data presented 
in lKim et al.l ( 19971 ) imply a mean separation of Lyman-a 



2 Note that 8 sub-DLAs (or "super-Lyman Limit Systems") 
were discarded so as to restrict the sample to genuine DLAs, 
even though they have Ovi coverage and adequate S/N. 
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forest lines (with logiV Hl = 13.1 - 14.3) of 560 kms" 1 
at 2=2.17-2.45, and of 360 kins" 1 at 2=2.71-2.30. We 
adopt an average separation of 400 kms -1 for the pur- 
poses of this calculation, ignoring the effects of clustering. 
The probability of a forest line falling in a ±20 krns" 1 
interval around z a b s is -Pbicnd ~ 40/400 ~ 0.1. The prob- 
ability of two blending lines falling at the same velocity, 
one in O vi A1031 and the other in A1037, or equivalently 
two blends falling in one O vi line at the velocity of two 
known components, is -Friend ~ 0.01, and so on for multi- 
ple lines. We find several DLAs with multiple O VI compo- 
nents identified; there is a negligible probability that they 
are all blends. In other cases there are only one or two 
components identified; these detections are less secure. To 
further justify our identification of Ovi in each DLA, we 
present detailed notes on each system in the appendix, 
which should be consulted together with Figure 1, show- 
ing the spectra. 

2.3. The resulting sample 

After applying our search criteria we are left with 12 O VI 
detections, of which three are at di>=c|,Zq SO — z a bs|/(l + 
^abs) < 5 000 krns" 1 from the QSO redshift, and so may 
be "associated" by the traditional definition (though see 
Sect. 5.1). C iv and Si IV are seen in every case where data 
is available, and Nv is detected in 3/9 systems with data. 
The Ovi detections cover the redshift range 2.07-3.08, 
median redshift 2.62, and are observed at S/N levels (per 
resolution element) in the range 23-41 at the Ovi wave- 
length. The S/N is typically over twice as high near Civ 
and Si iv. The DLAs span values of logA'ni from 20.35 
to 21.75, with metallicities [Z/ITjNeut (measured from the 
singly-ionized species in the neutral phase) also span- 
ning two orders o f magnitude between —2.59 and —0.49 



(jLedoux et al.ll2006l l I. The three systems with z a b s 



are toward Q0528-250, Q0841+129, and Q2059-360. 



c qso 



2.4. Incidence of Ovi absorption in DLAs 

We detect O vi in 12 of 35 DLAs (34%) with O vi coverage. 
There are no clear O vi non-detections, so the percentage 
of DLAs with Ovi could be 100%. However, we caution 
that O vi non-detections are difficult to make, due to the 
high density of forest lines; a conservative estimate of the 
fraction of DLAs with O vi is thus >34%. Analyzing only 
the intervening DLAs at >5 000 kins -1 from the quasar 
redshift, we find that 9 of 31 show Ovi, so the incidence 
of Ovi in intervening DLAs is >29%, whereas 3/4 DLAs 
at <5 000 krns -1 from the quasar redshift show Ovi. 



2.5. Measurement of absorption 

For each DLA we determined by eye the velocity range 
of absorption (labelled t>_ to v + ) in Ovi AA1031, 1037, 



Si IV AA1393,1402, Civ AA1548, 1550, Nv AA1238, 1242, 
and the neutral gaifl The exact range of ions available 
changes from system to system, due to blending and 
data availabil ity. Our neutral-pha se line for each DLA is 
adopted from lLedoux et al" ( 2006f ). who chose a line with 
0.1 < F(v)/ F c (v) < 0.6 in the strongest component to 
minimize saturation, where F c (v) and F(v) are the con- 
tinuum and observed fluxes as a function of velocity. 

We fit continua locally to all the lines of interest us- 
ing low-order (often linear) Legendre polynomials to re- 
gions of continuum free from absorption. We then in- 
tegrated the absorption in each lin e using the appar- 
ent op tical depth (AOD) method of ISavage fc Sembach 
(|199ll ). which returns measurements of the column den- 
sity N a = [3.768 x 10 14 /(A /)] T a (v)dv krns" 1 , where 
r a {v) = ln[F c (v)/F(v)], Ao is the rest wavelength in 
Angstroms, and ,f is the oscillator strength of the line 
(taken from iMortonl l2003h . The AOD method also re- 
turns measurements of the central velocity v of the ab- 
sorption using the first moment of the optical depth pro- 
file, v = f * VT a (y)dv/ /J + T a (v)dv. The treatment of er- 
rors arising from statistical (photon noise) and systematic 
(continuum placem e nt) ef fects is handled as described in 
Sembach fc Savage! (|l992l ). The UVES instrumental res- 
olution of 6.6 krns -1 (FWHM) ensures that lines with 
b > 4 krns -1 will be resolved, which is required for the 
AOD method to give accurate results. 

For each line we calculate Av, the velocity range con- 
taining the centra l 90% of the total integra ted optical 
depth in the line ( Prochaska fc Wolfei 1997b ). For cases 
where no absorption is seen at the DLA redshift in Nv, 
we present a 3u upper limit to the column density, using 
the 3cr upper limit to the equivalent width, and assum- 
ing any absorption to be optically thin. For cases where 
the line profile is saturated (which we define as where 
F(v) — <J F ( v j < in any pixel within the line profile), 
we present a lower limit to the column density, and an 
upper limit to the line width. 

We also used the VPFIT software packag^f] to derive 
the absorption line parameters for each high-ion absorber 
with Voigt profile fits. This provides an independent check 
of our absorption line measurements. The number of com- 
ponents fit is somewhat subjective; we use the minimum 
number necessary. We present the measurements using 
both the AOD and VPFIT techniques in Table 1. Note 
that the AOD method, which can be fully automated for 
high efficiency, and has the advantage of not requiring 
any assumptions about the component structure, gener- 
ally reproduces the same column density results (within 
the errors) as the VPFIT method. This indicates that the 
choice of u_ and w + by eye is not a significant source of 
error. 



Throughout this paper we use the convention 
[Z/H] =log(iVz /Nb )-log(JV z /JV H ) q • 



4 Note that the "neutral" phase is traced here through ab- 
sorption lines from singly-ionized metals; it is so called because 
hydrogen is predominantly neutral in this phase. 

5 Available at http://www.ast.cam.ac.uk/~rfc/vpfit.html 
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All velocities (in Table 1 and throughout this paper) 
are quoted relative to the redshift of the neutral gas in the 
DLA, which is measured using the central velocity of the 
stronge st component of the low-ionization metal absorp- 
tion in iLedoux et al.l (|2006l ). This provides a convenient 
rest-frame in which to discuss the kinematics of the ab- 
sorption line profiles. 

3. Results 

In Fig. 1 we show the absorption line profiles of all 
the DLA systems with Ovi detections. Fig 2. shows the 
New/ No vi ratios versus velocity for the ten DLAs in our 
sample with both O vi and C iv data. 

There is considerable variation in the appearance of 
the highly ionized absorption lines in the DLAs. The O VI 
absorbers range from cases with a single, optically thin 
component (e.g. toward Q2059-360 or Q2206-199) to cases 
with a series of saturated components (e.g. toward Q0450- 
131 or Q2243-605). The C IV profiles range from cases with 
one or two components spanning <100 kms -1 to cases 
with J>15 components spanning several hundred kms^ 1 . 
By no means is it clear that we are looking at a homoge- 
neous sample of objects. 

3.1. Ovi column density 

The logarithmic integrated O VI column density in our 12 
detections ranges from 13.66 to >15.15, with a median 
value of 14.77, a mean of 14.54, and a standard deviation 
of 0.59. The mean column density is similar t o that seen in 



the " thick disk" of O vi in the Milky Way (|Savage et al 
20031 ). where logAWi ranges from 13.85 to 14.78, with 
a median value of 14.38. A typical complete sight line 
through the Mi lky Way halo also shows O VI high- velo city 
clouds (HVCs; ISembach et al.l [20031 : iFox et alj|2006j) . so 
a better analogy for DLAs is the total Ovi column den- 
sity through the Milky Way halo over all velocities, which 
is typically 14.80. The average value of A^ovi in DLAs 
is slightly larger than both t he Large Magella nic Cloud 
average of (logAT vi)=14.37 (|Howk et al.ll2002l) . and the 
amount log ^0^=14.3 seen in starburst galaxy NGC 1705 



( Heckman et al.l 120011 ) . The values of JVhi in the neutral 
gas in the DLAs, between 20.35 to 21.75, compare to 
values in the Galactic IS M of between 19.6 and >22.4 
l|Dickev fc Lockman|[l990l) . with over two orders of mag- 
nitude of dispersion in each case. The similarity in the 
amount of O vi seen is intriguing considering that the 
metallicities in these DLA systems are typically forty 
times lower than the (~solar) value in the Galactic ISM. 
This implies that the total column density of ionized gas 
in DLAs is very large (see Sect. 5). 
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Fig. 2. Nctv/Nqvi ratio versus velocity. Each data point 
on this plot is a velocity pixel where both C IV and O vi 
are detected and unsaturated. jVciv is derived from the 
optical depth profile of A1548 unless this line is saturated, 
when we use A1550. iVovi is derived from A1031 unless 
this line is blended, when we use A1037. In nine out of ten 
cases (the exception is the DLA toward Q1337+113), this 
ratio is non-linear with velocity, showing that the O vi and 
C iv ion are generally not co-spatial. 



3.2. Comparison of high ion line profiles 

Throughout this Section we refer to Table 2, comparing 
the integrated kinematic measurements of the neutral- 
phase, C iv, and O vi absorption line profiles, Fig. 3, show- 



ing the column density, line width, and central velocity 
distributions for the components of Siiv, Civ, and Ovi 
in the DLA sample, and Fig. 4, comparing many measured 
absorber properties using scatter plots. 
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Fig. 1. VLT/UVES absorption line spectra of all DLA systems with detections of Ovi absorption. The tracer of 
the neutral gas is shown in the top panel, with the other panels showing all available high-ion data. In each DLA 
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Fig. 1. (-continued). All DLA systems with Ovi absorption. The lower three systems are at <5000 kms 1 from 
QSO. 
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Fig. 4. Comparisons between measured properties of O VI DLAs. TV is the total column density over all components, 
in units of cm -2 ; Av is in kms -1 . When a positive correlation between the two variables plotted is found (treating 
the upper/lower limits as data points), we display the Kendall rank correlation coefficient r and its significance on 
the panel, together with a least-square linear fit to the data (solid line). Dotted lines show one-to-one relationships. 
Saturated absorbers are shown with lower limits to log TV. The three DLAs at <5 000 kms -1 from the quasar redshift 
are highlighted in circles. 



The properties of the Siiv and Civ absorbers are 
strongly correlated. This is revealed by visual examination 
of the data, by the close similarity in their distributions of 
individual component line width and central velocities in 
Fig. 3, and by the strong correlation in their total column 
densities (Fig. 4, panel a). In general the Si IV and Civ 
thus arise in the same gaseous phase. 

There is a >3<r correlation between A^qvi &ud Nq w 
(Fig. 4, panel b). The total line widths of the two ions 
are also correlated (Fig. 4, panel c). Despite these correla- 
tions in the integrated properties, we do observe significant 
differences in the C iv and O vi absorbers at the compo- 
nent level. The Civ/Ovi column density ratios are non- 
linear with velocity in nine out of ten cases where data on 
both ions exist, often showing multiple peaks and troughs 
within the velocity width of each system (Fig. 2). This 
can be attributed to the O vi profiles being smoother than 
those of C iv, with fewer O vi components overall and no 
narrow O VI components observed. These non-linear ratios 
imply the Civ and Ovi arise in separate phases. Further 
evidence for O vi and C IV tracing different gas phases in 
DLAs is given by the histograms of the line widths of the 
individual components (Fig. 3, panel b), where strong dif- 
ferences are seen between C IV and O VI. Many narrow 
(b 55 10 kms -1 ) Civ components are observed, but most 
O VI components have b- values in the range 25-40 kms -1 . 



We report three DLAs with detections of both Ovi 
and Nv. One of these DLAs is close to the QSO redshift 
(toward Q0528-250) leaving two intervening cases where 
Nv is seen (toward Q0450-131 and Q2243-605). The Ovi 
column densities are highest (in fact, saturated) in the 
cases where Nv is present (Fig. 4, panel d). Profile com- 
parisons are hence not useful in revealing the O vi-to-N v 
relationship. 

3.3. High ions versus neutral species 



In Fig. 5 we compare the velocity range w_ to i>+ between 
O vi, C iv, and the neutral lines in each DLA, to display 
the larger range of velocity space that the ionized lines in- 
variably occupy. The mean total width (and its standard 
deviation) of all the Civ absorbers (including the lim- 
its from the saturated cases), (Avciv)=205±I25 kms -1 , 
is twice as broad as the average width of the neutral 
lines, (Au neu t)=102±88 kms -1 . Thes e results agree with 
the corresponding numbers reported in I Wolfe et alj (|2005l ) 
over a sample of over 70 DLAs: (Ai>c rv }=209±114 kms -1 
and (Au ncut )=II4±84 kms -1 . 

Individual C iv components are found at velocities of 
up to 400 kms -1 from the neutral gas (see the distribution 
in Fig. 3, panel c). We also observe offsets in the overall 
mean velocity v of high-ion absorption in each DLA. The 
mean and standard deviation of wciv is 65±78 kms -1 . 
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Fig. 3. Normalized histograms of the properties of high- 
ion components that comprise the DLA absorbers, as mea- 
sured using VPFIT. The number of components in each 
sample is indicated on the plot. Panel (a): Column density 
distributions. Panel (b): Line width distributions. Note 
that 6ovi > ^civ > ^Siivj i-e- the average component width 
rises with ionization potential. Panel (c): Central velocity 
distributions. Note the large dispersion in individual com- 
ponent velocities, at up to 400 kms -1 . 



I 



v (km s"') 

Fig. 5. Velocity range of absorption v_ to v+ for O vi, 
C IV, and the neutral gas in each DLA, showing the larger 
velocity dispersion in the ionized phases. Tick marks in- 
dicate the mean velocity of the C iv and O vi absorption. 
Partial blending limits the velocity range where we can 
measure the O vi absorption in several cases, marked with 
arrows. The range v+ — w_ is always larger than Av, our 
primary measure of line width; it is used here to show the 
distribution of ionized gas around the DLA redshift. 



These offsets can be interpreted as the signature of galactic 
inflows and/or outflows. 

In Fig. 4 (panels e to h) we directly compare our new 
measured high- ion DLA properties with the Hi column 
densities and neutral-phase metallicities. We find that the 
column densities and line widths of the high ions are each 
uncorrelated with the column density of H I in the neutral 
gas. However, the plasma properties do depend on the 
metallicity of the neutral gas. This is revealed by the ten- 
tative detection (98% confidence) of a [Z/HJ^eut vs ^Ovi 
correlation (panel g), and the tentative detection (also at 
98% confidence) of a [Z/Hj^eut vs Avcw correlation (panel 
h). Note we use Avciv and not Auovi to illustrate the 
high-ion kinematics since C IV line widths are far better 
measured. All DLAs in the sample with [Z/H]iMcut > — 1.5 
have both logiVovi >14.77 (the sample median value) and 



Ave iv >250km s 1 . A similar correlation is found between 
[Z/Hjueut an d Nqiv (not shown). 

The correlation between the high-ion line width 
and metallicity in DLAs follows the observed corre- 
l ation between the low-ion line width and metallicity 
(|Ledoux et alj|2006t fWolfe fc Prochaskalll998l ). However, 
these correlations likely have very different explanations. 
The low-ion velocity dispersion is thought to be dominated 
by gravity, and so its correlation with metal licity may im 



ply a n underlying mass-metallicity relation (ILedoux et al 
2006). The high-ion line width may be correlated with 
metallicity since both high ions and metals can be pro- 
duced as a result of supernovae. Further observations of 
DLA high ions are needed to investigate the existence of 
this correlation, which in principle can be used to probe 
the origin of plasma in DLAs. 



Fox et al.: Hot DLA halos 



9 



4. Ionization 

To create the ions O VI, N v, C iv, and Si iv from the pre- 
ceding ionization stages requires 114, 78, 48, and 34 eV, 
respectively. In principle, these ionization energies can be 
supplied cither by electron collisions in hot plasma or by 
extreme-UV photons. Either ionization process is plausi- 
ble in the vicinity of DLAs, due to the local presence of 
type II supernovae, accreting intergalactic gas, and the 
ionizing extragalactic background. 



4.1. Model-independent insights 

Many of the C iv and Si iv components have narrow line 
widths, b Si 10 kms -1 implying photoionization is the ori- 
gin for these components. Collisional ionization is ruled 
out since &crv < 10 kms" 1 implies logT < 4.86, at 
which temperature l ittle C IV would be produce d by colli- 
sions with electrons ( Sutherland fc Dopitalll993h : a similar 
argument holds for Siiv. These narrow components are 
particularly noticeable in the DLAs toward Q2243-603, 
Q0112-306, and Q2138-444. There are also broader Civ 
components present in the data with b w 30 — 40 kms -1 , 
e.g. toward Q1337+113 and Q2059-360. However, no indi- 
vidual narrow photoionized O vi components are found; 
the narrowest Ovi b- value is 14 kms -1 , with the ma- 
jority of cases over 20 kms -1 . For comparison, a ther- 
mally broadened O vi component at a temperature of 
3 x 10 5 K, where Ovi peaks in abundance in collision- 
ally ionized plasma, has a line width &ovi = 18 kms -1 . 
The line widths of the N v components are generally over 
20 kms -1 , though the sample is small. We cannot infer 
conclusively from the observed line widths whether O VI 
and Nv are photoionized or collisionally ionized, because 
of the possibility of non-thermal broadening or unresolved 
components. However, the fact that we clearly see pho- 
toionized C iv and Si iv components but do not see these 
components in O VI suggests that the O VI arises in a 
separate, collisionally ionized phase. Analogous situations 
(with narrow components seen in C IV but no t in Ovi 



are o bserved in the low halo of the Mil ky Way (IFox et al 



2003) and in Lyman Limit Systems (jKirkman fc Tvtler 
1999 . 



4.2. Photoionization models 

We ran photoionization mode ls using the code C LOUDY 
(version C06.02, described in iFerland et al.l["l998l ). to in- 
vestigate the viability of photoionization as an origin 
mechanism for the observed high ions in the DLA at 
z a bs = 2.076 toward Q2206-199. This absorber was chosen 
as an illustrative case. Our methodology involved running 
a grid of models at different gas density and total hydrogen 
column density, until we found a solution that reproduced 
the high- ion column densities (see Fig. 6). 

The following assumptions were made: (1) the plasma 
in DLAs exists in a uniform, plane-parallel slab; (2) the 
plasma has the same metallicity as the neutral gas; (3) 



the plasma has solar relative abundance ratios; (4) the 
plas ma is exposed to the extr agalactic background taken 
from lHaardt fc Madau ( 19961 ). with a mean intensity at 
912 A of J912 = 10 -21 - 5 ergs cm -2 s -1 Hz -1 sr -1 at z = 2. 



The first result is that in this test case (the DLA 
toward Q2206-199), there is no single-phase solution to 
the Si iv, Civ, and Ovi observations. There are solutions 
just to the C IV and Si iv, or to the C iv and O vi data, 
which we discuss here. A model to the observed Nsi rv and 
New, assuming [Si/H] = [C/H], requires log Ah = 20.3 and 
log U ~ —2.0, where U is the ratio of ionizing photon den- 
sity to gas density, corresponding to lognu = —2.7, or an 
over-density 5 — p/p « 200. The Civ and Ovi data are 
reproduced, if we assume [C/H] = [0/H], by a model with 
log Ah = 20.6, ionization parameter logU = —1.0, gas 
density log (?iH/cm -3 )= — 3.7, and linear size I — 630 kpc. 
However, using a value for the Hubble constant at z = 
2.5 of #2.5=220 h 70 kms- 1 Mpc -1 that follows from 
a standard (0^=0.7, S1m=0.3) cosmology, the Hubble 
flow broadening 6(kms -1 ) = iJ2.5^(Mpc) over a path of 
630 kpc would create lines with widths of «140 kms -1 , 
several times larger than the observed C IV line width (see 
spectra in Fig. 1). In order to arrive at a physically reason- 
able size, one r equires [O/C1>0, as i s the case for metal- 
poor halo stars ( Akerman et al.ll2004l ). In this case a given 
^Vciv/^Ovi ratio can be reached at a lower ionization pa- 
rameter, corresponding to higher density and smaller size. 



If local galactic sources of radiation were present 
in addition to the extragalactic bac kground (e.g. 
Howk fc Sembachlll999HVladilo et al.ll200ll ). then (assum- 
ing the shape of the radiation field does not change) 
a given ionization parameter would be reproduced at 
higher density and smaller cloud size, and consequently 
the Hubble broadening would be smaller. However, one 
would not expect a strong flux of O yi-ionizing phot ons 
from galactic sources of radiation fe.g. lFox et al.ll2005l ), so 
though local sources of radiation could produce the C IV 
they could not easily explain the O VI. 



Nv is not detected in the Q2206-199 DLA, in agree- 
ment with the model prediction. For those DLAs where 
Nv is present, one cannot use the Nv/ClV ratios to de- 
rive the gas density since the rel ative nitrogen abun dance 
in DLAs can be highly non-solar ( Pettini et al. 20021) . and 
also because the Nv may be collisionally ionized. We also 
note that two of the three DLAs at <5 000 kms -1 from 
the QSO (toward Q0841+129 and Q2059-360) do not show 
Nv absorption, whereas Nv is expected in intr insic sys- 
tems photoionized by radiation from the QSO (jHamann 
19971 ). This suggests that the influence of photoionizing 
QSO radiation is not that strong in these two cases. 
Indeed, the ionization properties of the three z a b s ~ z qso 
DLAs are no different than those of the intervening DLAs, 
so they may be representative of the whole population. 
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Fig. 6. Investigation of whether photoionization by the extragalactic background can explain the high ions in the 
DLA at z a bs = 2.076 toward Q2206-199, assuming solar relative abundances. The high-ion column densities (N, left 
panel) and their ratios (N\/ N%, right panel) are plotted as a function of ionization parameter U, the ratio of ionizing 
photon density to gas density, as determined using a series of CLOUDY runs. The observations are shown as large 
symbols. We find the C IV and Si iv observations can be reproduced by a model with total hydrogen column density 
logA^H=20.3 and gas density log nn = —2.7. This solution is physically reasonable. The Civ and Ovi observations 
can be reproduced by a model with logiVH=20.6 and log tt-h = —3.7. However, this solution implies a cloud size of 
«630 kpc; we rule this out since the Hubble broadening of such an absorber would be broader than the observed C IV 
and Ovi component widths. [O/C]>0 is required to produce a solution for Civ and Ovi that is not unreasonably 
large; collisional ionization models are also possible for Ovi. The log^ axis only applies to the case logiVH=20.6. 



5. Total ionized column density 

We can calculate the total hydrogen column density in the 
hot (O vi-bearing) and warm (C iV-bearing) phases using 



TyHot 



N , 



N Q , 



(O/H) 



Ncut 



N 



Warm 
Hll 



/ovi(0/H)not /ovi(0/H)Ncut (0/H)n o t 

iVciv _ iVciv (C/H) Ncut 



/C IV (C/H) Warm /dv(C/H)Ncut (C/H)\Var 



(1) 

(2) 



where fo vi = No vl / No is the fraction of oxygen atoms 
in the five-times- ionized state, and /crv = Ncw/Nq- 
Assuming solar relative abundance ratios in both the neu- 
tral and plasma phases, we have [O/H] Hot = [Z/H] Ncut and 

[C/H]\Varm=[Z/H]Ncut- By using Zhot.) Z\Varm, and ZNeut 

as shorthands for the absolute metallicities in the hot, 
warm, and neutral phases, we can write 



ArHot 

JV Hn 



jvrWarm 
JV Hn 



No, 



^Ncut 



/ovilO[ z / H l^(O/H) z Hot ' 

_ Npiv ZNcut 

" / Clv lO[ Z / H W( C /H) Zwarm' 



(3) 
(4) 



To evaluate Eqns. (3) and (4) for ea c h DL A we 
take (O/H) = 1CT 3 - 34 from lAsplund et alJ (l2004h and 
(C/H) =1Q- 3 - 6 ^ from lAllende Prieto et al.1 (|2002h . and 



[Z/H]Ncut from iLedoux et al.1 (|2006l ). We cannot directly 



measure the metallicity of the plasma phases so we assume 
ZNeut/ZHot and Z Ncut /Z Warm to be equal to one in our 
calculations. The remaining unknowns are /ovi and few, 
which will depend on the ionization mechanism. The max- 
imum value fo vi reaches in collisionally ioniz ed plasma is 
0.2 at 10 5 - 45 K (jSutherland fc Dopital[l993l) . The maxi- 
mum for /civ in collisionally ionized plasma is 0.29 at 
10 5 00 K, and our photoionization models also show a max- 
imum /civ=0.3 in photoionized plasma. We can therefore 
estimate lower limits to the hot and warm hydrogen col- 
umn densities using the maximum values for fovx and 
few- 

The implied hot hydrogen column densities are com- 
pared with the H I column density in the neutral phase 
in Table 3 and Figure 7, where we find \ogN^°f ranges 
from >19.5 to >21.1, and log N^f m ranges from >19.4 to 
>20.9. These lower limits are typically on the same order 
as the H I column in the neutral gas, though we observe a 
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Fig. 7. Comparison of the total ionized and neutral col- 
umn densities in DLAs. The top panel directly compares 
the total ionized column density A^n in both the hot 
and warm phases to iVm. In the bottom panel we show 
the ionized-to-neutral ratio Ahh/Ahi vs iVm. The dashed 
line shows the median value of N§°^/Nni from our nine 
intervening DLAs. Data points from the three DLAs at 
<5 000 kms -1 from the QSO are highlighted in large open 
circles. 

considerable dispersion (over two orders of magnitude) in 
the value of Nan/Nni- Note that if some Civ co-existed 
with the O vi in the hot phase, as is possible in collisional 
ionization models, then the true hydrogen column in the 
warm phase would be lower than that calculated here. 
Therefore we do not sum the H n contents of the C iv and 
O vi phases, since this could constitute a double-counting 
problem. 

6. Contribution of DLA plasma to the cosmic 
baryon and metal density 

The contribution of H I in DLAs to the cosmic density can 
be calculated using 



fi Hl (DLA) = 



CPcrit JlO 203 



10- 



iW(JV H i)dJV H i 



(5) 



where /(AThi) is the Hi column density distribution func- 
tion. r2 Hl has b een calculated as ~ 1.0 x 10~ 3 , fai r ly flat 
with redshift dProchaska. Herbert-Fort, fc Wolfel 120051: 
Peroux et aD [200j, 12005J). In Eqn. (5) the lower integra- 
tion limit corresponds to the definition of a DLA, and 
the upper integration limi t corresponds to the b reak in 
the power law reported by IProchaska et al. I (|2005l) ; above 
log A^h i=21.5 /(iVm) becomes much steeper and the con- 
tribution of DLAs at these column densities to il Hl be- 
comes small. S ince f (H i) has a power law dependence with 
slope a ~ — 2 ( Prochaska et al J 20051) . it can be shown by 
integrating Eqn. (5) that equal bins of log Ah i contain 



equal contributions to Q, i.e. all DLAs contribute toward 
the mass density. The contribution of the hot phase of 
DLA to the closure density can be expressed in a similar 
way: 



fl (DLA) 



CPcrit 



10- 



JVhi/(JVhi)cUVhi. (6) 



The ratio A/hii/A^hi m ay depend on A^hi, but since we 
only have nine intervening DLAs, we do not have enough 
data to characterize this relationship. Therefore we pro- 
ceed by using our median values (N§°* /A^hi) > 0.4 and 
(N^f m /A^Hi) > 0.2 measured over the nine intervening 
DLAs with O vj^. Now the contribution from the ionized 
gas simply becomes 



fi (DLA) = (Ah h/ATh (DLA) 



(7) 



evaluating to ft HotHn (DLA) > 4 x 10" 
f} WarmHlI (DLA) > 2 x 10~ 4 . The contribution 



to be augmented by fi z 

j^Warm H II 



and 
of the 

metal s in DLAs ng'(DLA) = (Z/H)» Hl (DL A) = 4 x 
10 -7 ( Ferrara. Scannapieco, fc Bergeronll2005l ) then needs 
HotHlI (DLA) > 2 x 10" 7 and 
(DLA) > 1 x 10~ 7 . 
Note that the plasma contributions to the baryon and 
metal budgets would rise if the value of /o vi in the DLA 
plasma were <0.2. The total density of metals estimated 
by integrating the star formation history of the Universe 
up to z = 2 is f^ FH w 3 x 10~^ (usin g the star for- 
mation history from lBouwens et al.ll2004l and assuming 
a metal yield y=l/42). Stars in galaxies ap pear to con- 
tain k 20% of the total ( Bouche et al.l 12006). the contri- 
bution from the ISM in galaxies (Hi in D LAs) is 



and the IGM contains a furthe r w 5 - 25% (|Pettinill200i 
Berg eron fc Herbert-Fortll2005l ). The remaining metals (s=s 
50% of the total) are yet to be found, h ence the for- 
mulation of the "missing metal problem" (jPettinil [l99l 
Ferrara et alH2005l ). If fa vi in the DLA plasma were as 
low as 3 x 10~ 3 , which is the case for plasma in CIE 



at 10 6 K ([Sutherland fc Dopitalll99a iGnat fc Sternberg 
20061 ) , then the O vi-bearing plasma around DLAs would 



constitute fl^ otHn = 1.5 x 10~ 5 and hence solve the miss- 
ing metals problem. In this event the majority of oxy- 
gen atoms would be ionized up to Ovn and Oviii. The 
lines of these ions are in the X-ray, but unfortunately 
searches for high-redshift DLA absorption in Ovn and 
O vni to confirm the presence of 10 6 K plasma are be- 
yond the capabilities of current X-ray satellites. A scenario 
with O vi existing as a trace ionization state in million- 
degree DLA coronae is consistent with the broader O vi 



lines in our sample, since bo^ 



30 kms in a ther- 



mally broadened line at 10 6 K. However, this scenario is 
not required by the data; it is merely one possible ex- 
planation. Indeed, the notion that all the missing met- 
als are in hot DLA halos is implausible, since metals will 



Using the mean rather than the median raises the 
Afj^f/ATm ratio to ml. However we use the median since it 
is less sensitive to the effect of false O vi detections. 
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also be found in both the neutral and ionized phases of 
other classes of QAL system. This incl udes the plasma in 
sub-DLAs and Lyman Limit Systems (|Kirkman fc Tvtler 
19971 ll999HD'Odorico fc Petitiear]|200ll : |Prochaska et al' 
20061). which m ay als o probe galactic halos (see Figure 1 
in lMaller et al.l [2003h . and could easily contribute to the 
missing metals. It will be necessary to search for and char- 
acterize the O VI phase in QSO absorbers over all ranges of 
Nji i to fully measure the quantity of baryons and metals 
hidden in hot galactic halos at high redshift. 



7. Discussion 

A full discussion 



of 



Mo & Miralda-Escude 


Rauch et al. 


1997; 


Mailer et al. 


200§) is 



models of the ori gin 

in DLA s (e.g. |mT 

19961: iKatz et alT 
Wolfe fc Prochaskal 



of 



1996; 



1996; 



2000b 



is beyond the scope of this pa- 
per. However, we briefly mention two leading theories in 
light of the findings in our new data. 



7.1. Accretion 

Numerical simulations of cosmological st ructure formation 
(|Cen fc Ostrikerlll999l iDave et alj|200l[ ) predict that gas 
falling onto diffuse large-scale structures will be acceler- 
ated to supersonic speeds and become shock-heated to 
temperatures of 10 5 to 10 7 K, creati ng a phase known 
as the warm-hot IGM (WHIM)0 The IDave et~aT1 (j200ll ) 
WHIM evolution models predict <10% of the cosmic 
baryons are in the WHIM at z = 3 (at over-densities 
10 < p/p < 30), but 30-40% are at z = 0. Our data 
show that collisionally ionized plasma at T > 10 5 K ex- 
ists around DLAs at z ^ 2, and that the baryonic con- 
tent of the plasma is at least of similar order to that in 
the Hi phase. This plasma could be identified with the 
high-density portions of the WHIM, since DLAs presum- 
ably trace deep potential wells within dark matter ha- 
los, and hence represent likely sites for accretion to occur. 
However, the presence of metals in DLA halos implies that 
the gas has been processed by star formation; accretion 
would bring in pristine, low-metallicity material. 



7.2. Supernovae and galactic winds 

Star formation in DLAs will lead to supernovae and re- 
gions of hot, shock- heated plasma, which will become 
outflows if enough energy is injected. Galactic outflows 
are a topic of considerable interest, not least because of 
the metal polluti on and energy inject i on they bring about 
in the IGM (e.g. IVeilleux et alJliool lAguirre et alj|2001 , 
20051: lAracil et al.ll2004| )~ 



Our DLA observations show O VI and C IV components 
displaced fro m the neutral gas by up to 400 kms ; for 
comparison, lAguirre et al.l (|200lh find that 300 kms -1 



7 Note that the term "warm-hot" here is at odds with our 
use of the terms warm and hot throughout this paper. 



winds at z > 3 can account for the metallicity of the 
Lyman alpha forest. Our data also show evidence at the 
98% level for both a [Z/H]N CUt vs Nq vi correlation and 
a [Z/H]Ncut vs Ave iv correlation. These correlations are 
consistent with a feedback scenario where supernovae in 
DLA disks produce both metals and kinematically dis- 
turbed volumes of interstellar plasma. 



8. Summary 

We have performed a successful search for O VI and other 
high ions in a VLT/UVES sample of 35 DLAs with data 
covering Ovi at S/N>10, and studied the properties of 
the highly ionized phase and its relationship to the neutral 
gas. We now list the principal observational results. 

1. We detect 12 cases of DLA Ovi absorption, of which 
three are at <5 000 kms -1 from the QSO. No Ovi non- 
detections are found, but considering the difficulty in 
making non-detections, a conservative estimate of the 
incidence of Ovi absorption in DLAs is >34%. All 
the DLAs show C iv and Si iv when data is available, 
and Nv is seen in 3/9 DLAs with Nv coverage. The 
Ovi absorbers are characterized by logVovi between 
13.66 and > 15.25 with a median value of 14.77. These 
average column densities are similar to the amounts 
measured in the Galactic halo, LMC, and starburst 
galaxy NGC 1705, even though the DLA metallicities 
(measured in the neutral gas) are typically one fortieth 
of the solar values. 

2. The Si IV and Civ profiles are very similar, suggesting 
these two ions are generally co-spatial. Many of the in- 
dividual C iv and Si iv components are narrow enough 
(b < 10 kms -1 ) to imply photoionization as the origin 
mechanism. 

3. In nine of ten cases the Ncw/Noyi ratio is non- 
linear with velocity through the profile, implying that 
O vi and C iv are not fully co-spatial. No narrow, 
photoionized O vi components are seen in the data. 
Furthermore, CLOUDY photoionization models show 
no single-phase solution to the Siiv, Civ, and Ovi 
column densities. Thus, the observations suggest the 
plasma in DLAs contains two phases: a photoionized 
phase seen in Si iv and C iv, and a hot phase seen in 
O vi and N v that may also contain additional amounts 
of Civ. 

4. C iv and O vi components are seen centered up to 
400 kms -1 away from the neutral gas, indicating in- 
flows or outflows. The total width of the C IV absorbers 
is typically twice as large as the total width of the neu- 
tral absorbers. 

5. Our data show tentative evidence (98% significance) 
for correlations between [Z/H]Ncut and logVovi, and 
between [Z/H]Neut and Awcrv- This can be explained 
if the DLA plasma is produced as a direct result of 
supernova-generated outflows. The high ion properties 
are uncorrelated to Nm in the neutral gas. 
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6. Using the observed A*ci V and Nq v1 measurements in 
the DLAs from thi s paper, the metallic ities and A% : 
measurements from Ledoux et al. (2006), conservative 
assumptions for the O VI and C iv ionization fractions 
(/ovi < 0.2 and few < 0.3), and assuming that 
the plasma has the same metallicity as the neutral 
phase, we calculate the total hydrogen column density 
in the O Vi-bearing and C iv-bearing phases. We find 
log iVjj °* ranges from >19.5 to >21.1, and log7V^ rm 
ranges from >19.4 to >20.9. These column densities 
of ionized plasma are typically of the same order as 
logA%i in the neutral phase. Using the median value 
to iVjfj^,' /iVm, we find that the total baryonic content 
of the Ovi phase is HotHn (DLA) > 4 x 10~ 4 and 
the total metal content is ftg otHlI (DLA) > 2 x 10" 7 , 
i.e. >1% of the metal budget at z « 2. If the tem- 
perature in the O vi-bearing plasma is « 10 6 K and so 
/ovi <C 0.2, then the Ovi phase of DLAs can make a 
significant contribution to the missing metals problem. 
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Table 1. High Ion Column Densities in DLAs 
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a Alternate (SIMBAD-compatible) names are given in parentheses. 

b Redshift of quasar determined from position of Lyman-a emission. 

c v- and v+ are the lower and upper velocity limits of absorption in kms" 1 . 

d Column density measured using AOD method. Strong and weak refer to the two lines of the doublet. 
e Column density measured using VPFIT. The number in parentheses indicates the number of Voigt components 
in the model. 

' AOD measurement made over a limited velocity range, due to partial blending. TV(fit), determined using 
unblended data from both lines, may be higher. 
9 Line completely blended. 
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Table 2. DLA kinematics: neutral vs ClV vs Ovi 
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a Redshift of the system determined using the velocity centroid of the strongest component in the neutral line. 

b Av is the velocity range enclosing the central 90% of the total integrated optical depth in the line (in kms -1 ). 

c v is the optical-depth weighted mean velocity of the line profile, v = JJ + VT a (v)dv/ r a (v)dv (in kms -1 ). 

d Lines saturated, so Av could be lower. 

e Lines partly blended and saturated, so no measurement is possible. 

* Lines partly blended, so v is unknown and only a lower limit to An is measured. 

9 No available C IV data. 

h Upper limits treated as data points in these calculations. 



Table 3. Total hydrogen column densities in Ovi phase, Civ phase, and neutral phase 
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>0.02 


Q1409+095 
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Q2059-360 


3.083 
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>19.9 
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>0.04 


Q2138-444 


2.852 
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Q2206-199 
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Q2243-605 


2.331 
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>14.98 
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" [Z/H]Nout 


and log 


ATki fromlLedoux et all (|2006l 


). Element Z 


is Zn if Zn n 


is detected, otherwise S or 


Si. 







b ^hi = ^Vo\i/[/ovi(0/H)Hot]; we present a lower limit since fom < 0.2. We assume [0/H]Hot=[Z/H]Neut. 
c N^n™ = ^civ/[/civ(C/H)warm]; we present a lower limit since few < 0.3. We assume [C/H] warm = [Z/H] Ncut . 
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Appendix A: Notes on individual systems 



Q0027-186 z, 



qso 



2.55 



abs = 2.402 

Ovi A1031 is blended below —50 kms -1 and above 
200 kms , but shows strong absorption components near 
0, 50, 125, and 175 kms -1 , as does the Civ profile. 
O vi A1037 appears to be unblended in the ranges —100 to 
60 kms' 1 and 180 to 220 kms' 1 , where the profile shows 
similar components to those seen in C iv. In our analysis 
we adopt log Vovi=15.09±0.03, based on the results from 
a Voigt fit to the unblended parts of each O vi profile. The 
AOD integration of Si iv A1393 yields a slightly higher col- 
umn density (by 0.11 dex) than A1402, so A1393 may be 
partly blended. 

QO 11 2-306 z qso = 2.99 z abs = 2.702 
Both Ovi A1031 and A1037 are blended at d < 
-100 kms' 1 . At v > -100 kms' 1 both Ovi lines show 
the same profile, which is clearly smoother than the C iv. 
We only present a lower limit to No vl in this system using 
a Voigt fit to the data at v > —100 kms' 1 . At least eight 
components are seen in the C IV and Si iv data. 

Q0450-131 z qso = 2.25 z abs = 2.067 
Both N v A1238 and A1242 show three components. Strong 
corresponding O vi is seen in A1031 and A1037, which have 
consistent optical depth profiles. C iv and Si iv data do not 
exist. The QSO is 17900 kms' 1 away from the DLA. Full 
details of the neutral gas i n this system are published in 



Dessauges-Zavadskv et al.1 ((2006) . 



Q0528-250 z qso = 2.77 z abs = 2.811 
This DLA is 3 600 kms' 1 more redshifted than the QSO 
and may be a QSO intrinsic absorber. Though the Ovi 
continuum placement is difficult due to a high density of 
lines, absorption features appear in A1031 at —150, —50, 
and 25 kms' 1 , the velocities of strong Siiv and Civ ab- 
sorption components. Our value for Novi derives from a 
four component fit to the data. The data show two very 
good detections of Nv with identical column densities 
measured in each line. 

Q0841+129 z qso = 2.50 z abs = 2.476 
This system is separated from the host quasar by 
2050 kms -1 . Though there are several blends near Ovi 
A1031 and A1037, two components appears in both lines 
with the same optical depth profile in the range —100 to 
kms' 1 . with the stronger component at —65 kms -1 . 
No Civ or Siiv data are currently available, and Nv is 
not detected. The neutral gas in this DLA is studied in 



Dessauges-Zavadskv et al.l (|2006l ) 



Q0913+072 z qso = 2.78 z abs = 2.618 
There are several blends near Ovi A1031 and A1037, but 
we verify Ovi absorption in the velocity range —22 to 
80 kms -1 since the A1037 profile shows the same column 
density profile as A 1031 in this interval. The C IV and Si iv 
data show at least four components over an interval of 
>200 kms' 1 . 

Q1337+113 z qso = 2.92 z abs = 2.796 
Ovi A1037 is blended, but the A1031 profile shows a sin- 
gle component with a similar profile to Civ, giving a 
Nciv/Nqvi ratio that is, within the errors, flat with ve- 



locity over a 60 kms 1 range. This is the only DLA in our 
sample where this ratio is flat. No Nv is seen. 

Q1409+095 z qso = 2.85 z abs = 2.456 
Ovi A1031 is blended, but the A1037 profile shows the 
same two components (central velocities of —19 and 
35 kms -1 ) as are seen in Civ and Siiv. In Civ and 
Ovi, the 35 kms -1 component is stronger than the 
— 19 kms -1 component. In Siiv and Sin, the —19 kms -1 
and 35 kms -1 components are of equal strength, showing 
that the Siiv follows the neutral gas in this case. Nv is 
not detected. 

Q2059-360 z qso = 3.09 z abs = 3.083 
This DLA is only 500 kms -1 away from the QSO. Ovi 
A1031 is blended; Ovi A1037 shows a clean detection of 
a single component, with a similar width to the Civ ab- 
sorption. Three components appear in the higher S /N C IV 
data, so the Civ/Ovi ratio is non-linear with velocity. 
This represents our highest redshift DLA O VI detection. 

Q2138-444 z qso = 3.17 z abs = 2.852 
O vi A1031 shows absorption between —300 and kms -1 , 
but is blended above kms -1 . The Ovi A1037 profile is 
identical to A1031 in the range —160 to 20 kms -1 , con- 
firming the detection of O VI, though the continuum place- 
ment is difficult. The C iv profiles broadly follows the O vi 
A1031 profile in the range —300 to kms -1 , although 
more subcomponents appear in the C IV data. 

Q2206-199 z q so = 2.56 z abs = 2.076 
This sight line shows a single O vi absorption compo- 
nent in both A1031 and A1037. Two components are 
revealed in the higher S/N Siiv and Civ data, sepa- 
rated by 20 kms -1 . The Ovi appears to be aligned with 
the weaker Civ componen t centered at 13 kms -1 . See 
Prochaska fc Wolfd ()1997al ) for a detailed analysis of this 
DLA; our total Civ and Siiv column densities arc in 
agreement with the earlier published values. 

Q2243-605 z qso = 3.01 z &hs = 2.331 
O vi A1031 and A1037 are each blended at different ve- 
locities but Ovi components centered near —350, —200, 
and —100 kms -1 correspond to features seen in Civ. A 
Lyman-7 line at z = 2.531 is identified at —270 kms -1 
in the O vi A1031 frame. We only present a lower limit to 
No vi in this system based on a 3-component Voigt profile 
fit to the unblended data. Nv is detected in both A1238 
and A1242, at velocities where the Civ and Ovi optical 
depths are highest. At least 14 components are seen in the 
Civ data over a velocity range of «500 kms -1 . Civ and 
Si iv follow each other in the series of absorption compo- 
nents in the range —225 to kms -1 , but the Civ pro- 
file shows many additional components in the range —500 
to —225 kms -1 that have no counterparts in Siiv. See 



Lopez et al.l ((2002) for a detailed analysis of this DLA. 
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